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Here we report our recent study on the spectral energy distribution (SED) of the high frequency BL Lac object 
Mrk421 in different luminosity states. We used a full-fledged x^-minimization procedure instead of more 
commonly used "eyeball" fit to model the observed flux of the source (from optical to very high energy), with 
a Synchrotron-Self-Compton (SSC) emission mechanism. Our study shows that the synchrotron power and 
peak frequency remain constant with varying source activity, and the magnetic field (B) decreases with the 
source activity while the break energy of electron spectrum (7;,^) a^nd the Doppler factor (5) increase. Since a 
lower magnetic field and higher density of electrons result in increased electron-photon scattering efficiency, the 
Compton power increases, so does the total emission. 



1. Introduction 

Active galactic nuclei (AGN) involve the most pow- 
erful, steady sources of luminosity in the Universe. It 
is believed that the center core of AGN consist of super 
massive black hole (SMBH) surrounded by an accre- 
tion disk. In some cases powerful collimated jets are 
found in AGN, perpendicular to the plane of accretion 
disk. The origin of jets are still unclear. AGNs whose 
jets are viewed at a small angle to its axis are called 
blazars. 

The overall (radio to 7-ray) spectral energy distri- 
bution (SED) of blazars shows two broad non-thermal 
continuum peaks. The low-energy peak is thought to 
arise from electron synchrotron emission. The lep- 
tonic model suggests that the second peak forms due 
to inverse Compton emission. This can be due to 
upscattering, by the same non-thermal population of 
electrons responsible for the synchrotron radiation, 
and s ynchrotron photons (S ynchrotron Self Compton: 
SScV lMaraschi et al.l (l992f . 

Blazars often show violent flux variability, that may 
or may not appear correlated in the different energy 
bands. Simultaneous observation are then crucial to 
understand the physics behind variability. 



2. x^-minimized SED fitting 

In this section we discuss the code that we have 
used to obtain an estimation of the characteristic pa- 
rameters of the SSC model. The SSC model assumes a 
spectrum for the accelerated electron density fc, which 
is a broken power law with exponents ni and 712. The 
minimum, maximum and break Lorentz factors for the 



DEF: SSC parameters initial values set-up 
calculate initial x'^ value, change parameters 
LOOP : 

calculate x for modified parameters 
if has increased: 

we are moving away from a minimum 
=^ change parameters, increase weight 
of steepest descent method and reset 
negligible decrease amount counter 
if has decreased: 

we are moving toward a minimum 

change parameters and increase 
weight of inverse Hessian method 
UNTIL: x' decreases by a negligible amount 



for the fourth time 



Table I The x minimization algorithm. 



electrons are usually called 7min. , 7max. and 7break re- 
spectively. The emitting region is considered to be a 
blob of radius R moving with Doppler factor 5 with 
respect to the observer in a magnetic field of inten- 
sity B. The model is thus characterized by nine free 
parameters. 

In the present work we have kept 7min. fixed and 
equal to unit, which is a satisfactory approximation 
already used in the literature. The determination of 
the remaining eight parameters has been performed 
by finding their best values and uncertainties from 
a minimization in which multi-frequency experi- 
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mental points have been fitted to the SSC spectrum 
modelled as in iTavecchio et al.l [l998j | . Minimization 
has bee n performed using the Levenberg-Marquardt 
method iPress et al. 1 |1994{ . which is an efficient stan- 
dard for non-linear least-squares minimization that 
smoothly interpolates between two different minimiza- 
tion approaches, namely the inverse Hessian method 
and the steepest descent method. For completeness, 
we briefly present the pseudo-code for the algorithm 
in table I. 

A crucial point in o ur implementation is that from 
ITavecchio et alj (l998j | we can only obtain a numeri- 
cal approximation to the SSC spectrum, in the form 
of a sampled SED. On the other hand, from table I, 
we understand that at each step the calculation of the 
requires the evaluation of the SED for all the ob- 
served frequencies. Although an observed point will 
likely not be one of th e sampled points coming from 
ITavecchio et al. I |1998I |. it will fall between two sam- 
pled points, so that interpolation can be used to ap- 
proximate the value of the SED^. 

At the same time, the Levenberg-Marquardt 
method requires the calculation of the partial deriva- 
tives of with respect to the SSC parameters. These 
derivatives have also been obtained numerically by 
evaluating the incremental ratio of the with respect 
to a sufficiently small, dynamically adjusted increment 
of each parameter. This method could have intro- 
duced a potential inefficiency in the computation, due 
to the recurrent need to evaluate the SED at many, 
slightly different points in parameter space, this being 
the most demanding operation in terms of CPU time. 
For this reason we set up the algorithm to m inimize 
the number of calls to ITavecchio et al.l 1998j | across 
different iterations. The x^ fit during different itera- 
tions are shown in Fig. 1. 



3. Application and results 

In order to study the behavior of parameters with 
source activity, we choose Mrk421 (table II), consid- 
ering the larger availability of MWL data sets and the 
lower redshift, hence less uncertainty after EBL cor- 
rection of VHE data. The fitted SEDs are shown 
in Fig. 2. 

In addition to the x^ test, we also checked the good- 
ness of the fit using the Kolmogorov-Smirnov (KS) 
test. Considering the occurrence of different physical 
processes (synchrotron and inverse Compton, at sub- 
stantially different energies), and the different quality 



^The sampl i ng o f the SED function coming from 



ITavecchio et al.1 Il998ll is dense enough, so that, with respect 
to other uncertainties, the one coming from this interpolation 
is negligible. 



Table II Data sets used in this study. The observation 
period of each state can be found at Fig. 2. 



State 



Instruments 



References 



1. XMM-Newton 
Whipple, MAGIC 

2. XMM-Newton 
Whipple, MAGIC 

3. KVA, WIYN, RXTE 
Whipple, HEGRA-CT 1 

4. Boltwood, RXTE 

Whipple 

5. Havard-Smithsonian 

RXTE, Whipple 

6. XMM-iVewton 

VERITAS 

7. Havard-Smithsonian 

RXTE, Whipple 

8. WEBT, Swift 
RXTE, VERITAS 

9. Boltwood, RXTE 

Whipple 



Acciari et al. [2009 



Acciari et al. [2009 



Rebillot et al. [2006 



Blazeiowski et al. [2005] 
Fossati et al. [2008] 



Acciari et al. [2009 



Fossati et al. [2008 



Donnarumma et al. [2009 



Blazeiowski et al. [2005 



of low- and high-energy data, we used a piecewise KS 
test, i.e. we applied the KS test separately to low- 
and high-energy data. Then the KS test always con- 
firms that the fit residuals are normal at 5% confidence 
level. 

Our results suggest that in Mkn421, B decreases 
with source activity whereas 7break and S increase 
(Fig. 3 top). This can be interpreted in a frame where 
the synchrotron power and peak frequency remain 
constant with varying source activity by decreasing 
magnetic field and increasing the number of low en- 
ergy electrons. This mechanism results in an increased 
electron-photon scattering efficiency and hence in an 
increased Compton power. Other emission parame- 
ters appear uncorrelated with source activity. In Fig. 3 
(bottom) , the B-7broak anti-correlation results from a 
roughly constant synchrotron peak frequency. The B- 
S correlation suggests that the Compton emission of 
Mkn421 is always in the Thomson limit. The S-jhicak 
correlation is an effect of the constant synchrotron 
and Compton frequencies of the radiation emitted by 
a plasma in bulk rclativistic motion towards the ob- 
server. 

More detailed description of this work is available 



More detail ed aescr i ption ^ 
in lMankuzhivil et al. I |2011 |. 
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Figure 2: Mrk421 MWL data sets and corresponding minimized SED fits. 
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Figure 3: Top. Variations of the fitted parameters - B, 5, and 7br - as a function of luminosity. The other SSC 
parameter s show a s catte r plot with the luminosity. Parameters and its uncertainty can be found at 
iMankuzhivil et al. I [201l| . Bottom. Correlations between B, 5, and 7br. 



References 

Acciari, V. A., et al. (VERITAS Collaboration) 2009, 

ApJ, 703, 169 
Blazcjowski M. et al. 2005, ApJ, 630, 130 
Donnarumma I. et al. 2009, ApJ, 691, L13 
Fossati, G., Buckley, J.H., Bond, I.H., et al. 2008, 

ApJ, 677, 906 
Mankuzhiyil, N., Ansoldi, S., Persic, M., Tavecchio, 



F. 2011, ApJ, 733, 14 
Maraschi, L., Ghisellini, G., & Celotti, A. 1992, ApJ, 
397, L5 

Press, W.H., et al. 1992, Numerical Recipes (Cam- 
bridge: Cambridge University Press) 

Rebillot, P. et al. 2006, ApJ, 641, 740 

Tavecchio, F., Maraschi, L., & GhiseUini, G. 1998, 
ApJ, 509, 608 (T98) 



eConf C110509 



